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Abstract: We give background material and some details of calculations for two 
recent papers [1,2] where we derived a path integral representation of the transi- 
tion element for supersymmetric and nonsupersymmetric nonlinear sigma models 
in one dimension (quantum mechanics). Our approach starts from a Hamiltonian 

H(x,p, ip, ij)*) with a priori operator ordering. By inserting a finite number of com- 
plete sets of x eigenstates, p eigenstates and fermionic coherent states, we obtain 
the discretized path integral and the discretized propagators and vertices in closed 
form. Taking the continuum limit we read off the Feynman rules and measure of the 
continuum theory which differ from those often assumed. In particular, mode regu- 
larization of the continuum theory is shown in an example to give incorrect results. 
As a consequence of time-slicing, the action and Feynman rules, although without 
any ambiguities, are necessarily noncovariant, but the final results are covariant if 
H is covariant. All our derivations are exact. Two loop calculations confirm our 
results. 



1 Introduction. 

The subject of path integrals in curved space is arcane, complicated and contro- 
versial [3]. In two recent articles we have considered one-dimensional (quantum 
mechanical) path integrals, and found an exact path integral representation for the 
transition element [1,2]. For the bosonic case it is defined by T(z, y; (3) = 
< z\ exp — &H(x,p)\y > where \y > and < z\ are position eigenstates. The classical 

1 to appear in proceedings of workshop on gauge theories, applied supersymmetry and quantum 
gravity, Leuven, Belgium, July 10-14, 1995, and strings'95, USC, March 13-18, 1995. 



Lagrangian is given by L d = ^gij(x)x l x^ , but both the quantum Hamiltonian and 
the action in the path integral deviate substantially from L d . We begin by assuming 
that H has a given a priori operator ordering. By inserting N — 1 complete sets of 
x eigenstates and N sets of p eigenstates, one finds a discretized phase space path 
integral, from which one can derive (as we shall indeed do) discretized propagators 
and vertices in closed form (by coupling to discretized external sources). In the 
continuum limit one obtains then a Euclidean path integral of the form 

J dpdxJ-f^- HM}dt (1) 

with well-defined propagators, vertices, and rules how to evaluate Feynman graphs. 
The last result is the most important: these rules are new and differ from what is 
usually assumed. 

If one would bypass a detailed analysis of the discretized case, one might expect 
that L = ipq — H(p, q) is covariant if H is a covariant operator. (For example, 
if H would commute with the supersymmetry generators, one might expect that 
after integrating out p, the resulting actions are the supersymmetric actions one 
encounters in the literature). This is incorrect: the action needs noncovariant terms 
of order h and h 2 (but not beyond) in order that T be covariant. 

Another source of puzzlement might be the observation that actions of the form 
L = ^gij{x)x l x^ contain double-derivative interactions, leading to linearly divergent 
graphs by power counting. On the other hand, it is well-known that quantum 
mechanics is a finite theory. It would seem strange (and is, in fact, incorrect) 
to require that normal-ordering removes divergences: where would normal ordering 
come from? The resolution of this paradox will be the presence of new ghosts, closely 
related to the factors g 1//2 5(0) which Lee and Yang found in a careful treatment of 
the deformed harmonic oscillator [4], and which we shall call for this reason "Lee- 
Yang ghosts" . (They were first introduced by Bastianelli, and in a more covariant 
form by him and one of us [5]). 

The propagator for the quantum deviations of a scalar q(r) (where r = t//3) on 
the interval (—1,0) with boundary conditions q(— 1) = q(0) = is proportional to 

A(<7, r) = A F (a - r) + (ar + ±r) (2) 

where A F (a — r) = — T)9(a — t) + \(t — a)9(r — a) is the translationally invariant 
solution of d 2 .A F (a — r) = 5(a — r) while the terms ot + |(X+|t enforce the boundary 
conditions. In Feynman graphs, contractions between x and x, and between x and 
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x, produce 9 (a — r) and 5(a — r), and the problem arises how to evaluate products 
of several 5 (a — r) and 9 (a — r). Mathematically, various consistent but different 
definitions of these products of distributions can be given [6], but physically (1) 
should have an unambiguous meaning, and the problem is to find the correct rules. 
Consider as an example 

J° J° 5(a - t)9{o - t)9(<j - r)dadr (3) 

One might expect that, since 5 (a — t) = d a 9(a — r), the result is equal to 

| J-i J-i d a [9(a — r) 3 ]dadt = |. However, the correct result is \ as the discretized 

approach shows. 

Another source of ambiguities are the equal-time contractions, in particular 
< x(t)x(t) >. Are these the limit of < x(ti)x(t 2 ) > for T\ — > r 2 , and what to 
do with the resulting 5(0)? In field theory, equal-time contractions are a priori 
undefined, and one needs a symmetry principle to fix them, but here everything has 
been specified from the beginning, so ambiguities in equal-time contractions are not 
allowed. 

Yet another worry would be the perennial headache called "the measure" . Since 
the path integrals correspond to a one-dimensional quantum field theory on a finite 
time segment, one might expect some factors in the measure to be present, like det g 
to some power at the end points. In fact, in Hamiltonian quantization of gravity 
in higher dimensions, such factors (and factors involving goo) are present. We shall 
see that also in our case there are nontrivial measure factors. These are usually 
omitted in calculations with nonlinear sigma models, but are crucial to obtain the 
correct results. 

What do we exactly mean by "correct results"? "Correct results" means for 
us: the results which agree with < z\exp— j-H\y >, no more and no less. This 
matrix element can be straightforwardly evaluated order by order in (3, without 
encountering any divergences or ambiguities, simply by expanding the exponent 
and inserting a complete set of p-states 

T(z,y;/3) = J < z\exp--H \ p >< p\y > d 4 p (4) 

Moving all x j to the left and all pj to the right, keeping track of commutators, one 
obtains c-number results. All terms with a given number of commutators, say s, 
are of a given order in j3 (see below), and although for each s an infinite number of 
terms contributes, one can sum the infinite series for fixed s in closed form. Thus 
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T(z, y; (3) is a finite and unambiguous Laurent series in f3. Our task is to find a path 
integral which reproduces these terms order by order in loops (in the path integral, 
(3 and h appear only in the combination j3h, so (3 also counts the number of loops on 
the worldline). One could, of course, reject this Hamiltonian starting point, and try 
to devise a more covariant way of defining path integrals in curved space which does 
not require a given Hamiltonian H(x,p). In particular, whereas in our approach 
we encounter noncovariant "midpoint rules" like Xk+i/2 — \{ x k+i + x k)i one might 
hope that covariant midpoint rules (the middle of a geodesic, for example) might 
lead to a completely covariant treatment. All we can say is that we have found a 
path integral which yields the correct results and which straightforwardly follows 
from T(z,y;j3) by inserting complete sets of states, whereas the more covariant 
approaches have so far not been able to reach the same results^. 

There are good practical reasons for starting with the Hamiltonian matrix el- 
ement < z\ exp — (3/TiH \y >, rather than a covariant configuration space starting 
point. When one calculates anomalies in n-dimensional quantum field theories, one 
can rewrite these anomalies as products of matrix elements of the Jacobian times 



The extension to include fermions is straightforward. So, quantum mechanics en- 
ters via T, and the reason one wants to rewrite T as a path integral is that the 
calculations of anomalies are much simpler in the path integral approach than in 
the Hamiltonian approach. We shall not discuss these applications to anomalies 
here, but refer to [2]. 

We shall also not discuss fermions in detail here (again, see [2]), but only say 
that one can introduce bras and kets \q > and < fj\ which are eigenstates of the 
fermionic annihilation and absorption operators ip a and xpl, respectively. These 
states are coherent states \q >= (exp^r] a )|0 > and < fj\ =< 0| expfj a ip a with rf 
and fj a independent Grassmann variables (so not related by complex conjugation, 

2 Covariant techniques to evaluate path integrals do exist, but can only be used when the path 
integral has additional special properties, for example if the semi-classical approximation is exact, 
or if the theory has additional symmetries so that one can use localization techniques. However, 
such techniques do not extend to arbitrary path integrals. In order to set up perturbation theory 
one has to choose a decomposition S = + S' nt , and in the case of a general a model it is not 
true that S^ ' and S lnt can be chosen to be separately covariant. (The kinetic operator in the full 
quadratic term is covariant, but not invertible in closed form[3].) Thus, Feynman rules based purely 
on are necessarily noncovariant. 



T[7], 
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only satisfying / dr/i] = J dfjfj = 1). This is the fermionic equivalent of the "holomor- 
phic" representation for bosonic systems [8]. We have found it simplest to use the 
x,p representation for coordinates, but for the fermionic part the iftjift^ representa- 
tion is by far the most natural. (One could, however, also construct a kind of x,p 
representation for the fermions) . The fermionic coherent states satisfy completeness 
relations and for N = 2 supersymmetric systems there is really no major obstacle 
to construct T(fj, z,rj,y; (3): one combines ift^a = 1,2) into ift a = (iftl +#2)/v / 2 and 
iftl = (ipi — iift a )/\/2. However, for iV = 1 supersymmetric systems, one has n Majo- 
rana fermions ift a (a = 1, . . . , n) satisfying the Dirac brackets {ift a , ift b } = 5 ab . To con- 
struct a vacuum and coherent states, we need rather ip A and ip\ satisfying ip A \Q >= 
0, {4) A ) il) B } = 0, {i>A,i>B} = and {ip A -,i^ B } = S A B - This can be achieved in two 
ways: 

(i) by fermion doubling, namely adding a second set of free fermions tpjj which do not 
appear in the Hamiltonian but which are used to construct ip A and ip\ as ip A = + 
iip i j) I V% and ift\ = (iftj — iift i j) I y/2. Here iftj denotes the original set of fermions ip a . 

(ii) by fermion halving, namely as follows 

^ = (^ A ~i + # 2A )/ v / 2, ift\ = (ift™- 1 - iift 2A )/V2 (6) 

The vacuum and Hilbert space are different in both cases, and one finds different 
results for the propagators and transition elements (!), but the anomalies come out 
the same. This is as expected: in traces over Hilbert space, differences created by 
choosing different vacua should cancel. 

In the next sections we shall give some details of calculations which are only 

5 



briefly summarized in [1,2]. A flow chart of the main ideas is given above. 



2 Weyl ordering and an extension of the Trotter 
formula. 

For definiteness one may consider a particular Hamiltonian H E which plays a role 
in the calculation of anomalies 

H E = \g ' 7',//' V'/',// ' ' (7) 

This operator is Einstein invariant, but we stress that our results hold for any other 
operator with two p's. Inserting (N — 1) complete sets of x eigenstates and N com- 



plete sets of p eigenstates, using the completeness relations / \x > y g(x) < x\d n x = 
J \p >< p\d n p = I, we find an expression for T(z, y; /3) in terms of N kernels 

T(x k ,Xk-i;e) = J <x k \exp-^H\p k ><p k \x k - 1 >d n pk, (8) 

where x N = z, x = y and e = (3/N. We rewrite if as a Weyl ordered operator. For 
a polynomial in p's and x's, the corresponding Weyl ordered operator is obtained 
by expanding (a^pj + PiX l ) N and retaining all terms with a particular combination 
of a's and /3's. It follows that (g lj PiPj)w — \PiPj9 lj + \Pi9^Pj + \g lj PiPj, and by 
evaluating H — {\g l ^PiPj)w one finds extra terms of order ft 2 

H = l(9 ij Pi Pj )w + y (T^T*^ + R) (9) 

(on a sphere R < 0). Other operators H in which the p's appear less symmetrically 
will in general also lead to extra terms of order h. For a polynomial one may prove 

(TTl \ 
x m ~ e p r x e , and it follows that 

m 1 / m \ f 

< z\{x m p r ) w \ y > = ^ I £ j J zm ~ e p r y e < Ap >< p\y > d n p 

r / z _|_ y\ m 

= J < z\p> f j p r < p\y > d n p (10) 

This shows why orderings like Weyl ordering are very convenient: one can replace 
a Weyl ordered operator by a function, simply by substituting p — > p, x — > |(z + y), 

and t/iis an exact result. 
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However, Weyl ordering and exponentiation do not commute, (exp —$H) W ^ 
exp-^(Hw), and whereas Hw was easy to write down, a closed expression for 
(exp— tH)w cannot be written down. One expects, however, that a suitable ap- 
proximation of the kernels, containing only terms of order e suffices. Here one 
stumbles upon a problem: it might seem that p is of order e -1 / 2 due to the term 
exp — |ep 2 in the action. Expansion of exp —%Hw would contain terms of the form 
e s p r f(x) for which s > 2 but which would still be of order e. We are now going to 
give an argument that p is of order unity, and therefore only the terms with one 
explicit e need be retained. Hence, we will use as kernel exp — ^Hw (\(xk+ Xk-i) , Pk) ■ 
In other words, the Trotter-like approximation 

< x\ exp — — H\p >~< x\l — —H\p >= (1 — —h) < x\p >=< x\p > exp —~h (11) 
n n n n 

is still correct, but h is not simply < x\H\p > as in the usual models with H = 
T{p) + V(x), but rather it equals Hw at the midpoints. To obtain this result, we 
note that the kernels are exactly equal to 



d n pke m-^-v* e-i* (x k „ 1/2 ,p k ); - —\ (12) 

V /W Xk-l/2 = 2 \ X k + X k -l) 

The difference between (exp —$H) W and exp — $(H W ) consists of two kinds of terms 

(i) terms without a p ; these are certainly of higher order in e and can be omitted 

(ii) terms with at least one p. 

The crucial observation [9] is now that the phase space propagators < Pk,iPi,j > and 
< Pk,i^k+i/2 > are both of order unity, and not of order e _1 and e -1 / 2 , respectively. 



A formal proof is given in equation (p5|) . However, already at this point one might 
note that the pp propagator is not only determined by the gpp term but also by 
ipAq. Completing squares, it is the p' = (p — iAq/e) which is of order e^ 1 / 2 . In the 
pp propagator the singularities of the p'p' and AqAq propagators cancel each other. 
As a consequence, the pp and pq propagators are of order one, and this proves the 
Trotter formula also for nonlinear sigma models. 

The final result is that one may use exp —&Hw(^(xk+Xk-i),Pk) as the kernels of 
the path integral. If one would not have used Weyl ordering, but simply computed 
/ < Xk\ exp — %H\pk >< Pk\Xk-i > d n pk keeping all terms of order e, one finds terms 
proportional to Rij(x~k-i/2)Ax k _ 1 / 2 Ax J k _ 1 / 2 where Rij is the Ricci tensor. These 
terms do not correspond to a local action (they are of the form / Rijif&dte). These 
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nonlocal kernels will yield the correct answer for the path integral, but Feynman 
rules for nonlocal theories are a headache, and it is not clear whether a truncation of 
these kernels to a local action exists which yields the correct answer. Weyl ordering, 
on the other hand, does lead to local kernels which are very easy to construct and 
which yield the correct result. This demonstrates the usefulness of Weyl ordering. 



3 Discretized propagators and new Feynman rules. 

If we keep the p k j and x\, as integration variables we obtain a discretized phase space 
path integral, but if we integrate the p k j out, we get a discretized configuration 
space path integral with N factors [det • In the continuum limit Lee 
and Yang wrote these determinants as exp |5(0)£r In gijdt and treated the exponent 
as a new term in the action [4]. We first discuss these configuration space path 
integrals. 

For the same calculational advantages as in the case of Faddeev-Popov ghosts 
in gauge theories, we exponentiate these determinants by ghosts, but whereas the 
Faddeev-Popov determinant needs a pair of anticommuting real ghosts/antighosts, 
here the square root of the determinant requires an extra third commuting real 
ghosts. One obtains then the following path integral [5] 

N-l 

n ( det 9ij(x k+ i /2 )) i/2 = 

k=0 

a j db k+1/ 2dc k+1/2 da k+1/2 exp-g ij (x k+1/ 2) [&fc+i/ 2 4+i/2 + a k+i/2 a { + i/ 2 ] ( 13 ) 

where a is a constant which can easily be determined. We decompose the x k into a 
sum of background parts and quantum parts, x\ = x{ g k + q{, and the action S into 
a free part 



iV — 1 

^ (0) = E o9a(z) (Aql +1/2 Aql +1/2 + V k+1/2 c? k+1/2 + ai +1/2 a{ +1/2 ) (14) 

fc=0 z 

and an interaction part S mt (the rest), requiring that the background fields satisfy 
the (discretized) equation of motion of and the boundary conditions. Hence 
Qn = % = 0. With the continuum limit in mind we parametrize the discrete q k using 
continuum modes of 

^=£^4sin^;fc = l,...,7V-l (15) 



m=l 



The Jacobian for x k — > q k — > r m is unity. We then couple Aq 3 k+1 ^ 2 and q J k+1 / 2 to 
discretized external sources 



JV-1 



^(sources) = - ^ { F k+i/2,Aq 3 k+l/2 + eG fe+1/2 , g£ +1/2 ). (16) 

fc=0 

Similarly we introduce discretized sources for the ghosts b, c and a. We then com- 
plete squares and integrate over the discrete variables r m J , b 3 \+i/2, c 3 k+i/2i aJ k+i/2- 
The result is a functional quadratic in external sources which will yield the dis- 
cretized propagators. We first quote the result and then give details of the calcula- 
tion. By differentiating twice w.r.t G one finds 



< ^+1/2^+1/2 >= 



-(* + \)(* + + (£ + \)0 k/ + (k + 1/2)9^ 



;i7) 



where 9 k ^ is the discretized 9 function (9 k ^ = if k < £, 9 k ^ — | if k — £ and 9 k ^ = 1 
if k > £). In the continuum limit this becomes 

< ?' > = -/3hg*(z)A(a, r) ; -1< a = ^±4 - 1< 

A(cr,r) = <t(t + 1)0(<7 - t) + t(<t + 1)0(t - <r) (18) 

Similarly 



< g i fe+ i/ 2 Ag^ +1/2 > = e%^'(z) 



2 +9 k 



N 



<Aq\ +1/2 Aq 3 t+l/2 > = € ^(z)[-l + N5 k , e ] 

2 1 

< V 11+1/20*1+1/2 > = --ftg l3 {z)5k,t\ < a l k+i/2a ] i+1/2 >= -hg v {z)5 k/ (19) 

These results show that in the continuum limit 9(<j — t) = 1/2 at o = r, 5(<j — t) 
is a Kronecker delta even in the continuum theory and not a Dirac delta, and they 
define equal-time contractions. For example 

< 9* (o-)q j {<r) > + < b l (a)c j (a) > + < a\a)a 3 {a) >= -(3hg i3 {z) (20) 

We see how the ghosts remove divergences, but we also see that a well-defined finite 
part is left which in a less rigorous approach might have been missed. Terms in 
Feynman graphs with more than one S(a — r) are eliminated by the Lee- Yang ghosts 
whereas products of one 5(a — r) and any number of 9(a — r) are evaluated by still 
defining 5 (a — r) in the continuum case to be a Kronecker delta. 

9 



4 Derivation of the discretized propagators. 



The orthonormality of the matrix O k m = ^ sin ^kmir in flT5| ) follows from the 
formula 2 sin a sin (3 = cos(a — (3) — cos(a + j3), and for — 2N < p < 2N 



N-l 

E 

m=l 



COS : 



prrnr 1 
~N~ ~ 2 



N 



1 1 



E e4pWAr - 1 - R p I = ns p,o - § - §(-)" !21! 



■ m=-iV+l 



The free action S^ ** is diagonal in r's since the O's in Aq k _iAq k _i in (|ll]) appear 



as 



A' 



iV 



(O k m - O fc _D (O fc n - O k ^ n ) = 25 mn - Ok m (O k ^ n + O k+l n ) (22) 



k=l 



k=l 



and using sin a + sin (3 = 2 sin |(a + /3) cos ~(a — /?), the orthogonality of O k m leads 
to 

^^Lww ( 1 - cos ^ [ ) 



m=l 



(23) 



Adding 5" (sources) to 1 we find 



n N-l 



Z[{F},{G}] = I 1] II ^-exp-- |^ + 



j=l m=l 



7V-1 n 



E E { -( F k-l/2,j ~ Fk+l/2,j) + g ( G k-l/2,j + G k+ l /2 ,j) > < E V "at 



fc=l j=l 

Completing squares and integrating over r m J yields 
w-i (7r e ^)»/2 



f JV-1 



r m J sin 



n 



i det^(^) 1 /2(i _ cos^vr)™ 



/2 



exp 



■jv-i 



4(1 - 



I m=l 



krmr 
~N~ 



i 4(1- cos ^) 



^(F,G) S 



2/2 m7r 1 m7r 

^(F,G) = iV Sm 27V cos ( fc + ^>W Fk+1 ^ 

12 mn 1 ^ 1 m7r 

+ W COS > Sin feH ) Crfc+1/2 7 

V A 2 A ^ V 2 A + 7 ,J 
The square of Q is, of course, taken with g^(z). 



(25) 



The easiest propagator to compute is < qq >. By differentiation w.r.t. -F fe+1 / 2 ,i 
and Fi + i/2,j one finds that the square of sin ^ cancels the factor (1 — cos ! ^ L ) in the 
denominator, and using cos a cos (3 = | cos(a + /3) + | cos(a — /3), one must evaluate 
the sums J2m=i °f cos (^ + ^ + l)mn/N and cos(/c — £)mir/N, for which one may use 
(FH). The result is given in (^9|). 
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Next we consider the < qq > propagator. Differentiation w.r.t. Gk+i/2,i and 
Fe+1/2,3 leads to a product cos — cos(£ + — sin (A; + — sin The last factor 
partly cancels the denominator 1 — cos One must then evaluate the series 

1 AT-1 

_ ^ ^ ((^ m "I - C~ m )(C^' Cm "I - C^^ _ "I - (^~ 2 ^ m ) _|_ £~ (2£+l)m\ ^26) 

4 m=l 

where we defined ( = exp J^. We write this series as a sum of four series, and 
combine terms pairwise such that we can use ([21]). The first two series start with 

£(2fc+2£+2) m and £(2fc+2£)m respective ly 5 an d run t iH £(-2fc+2^+2)m and ^{-2k+2i) m ? while 

the last two series we write in ascending order such that they start with ^(- 2k - 2e ~ 2 ) m 
and Q(- 2k - 2i ) m anc i run till ^( 2k -2£-2)m an ^ ^(2k-2i)m^ reS p ec tively. The terms in the 
first and third series are pairwise combined using (^lj), and similarly the terms in 
the second and fourth series. One finds then 

k-\-t-\-\ -i k+l 

7 E (-l-H P + 2iV5 Pi0 ) + - J! (-l-H P + 2iV5 p , ) (27) 

The terms with (— ) p cancel. In the remainder one may distinguish the cases k > 
£,k < £ and k — I. One finds then \ [-2{2k + 1) + 2NS k >i + 2N5 k>i \. So the result 
is proportional to (k + 1/2) /N + \5k>t + \$k>£ which agrees with (|19|). 

Finally we consider the qq propagator. This is the most complicated one. Dif- 
ferentiation with and Ge + i/ 2 j leads to the series 

^ run / sin(fc + l/2)^ \ / sin(l + i)- \ 

m ^ 1 (C ° S 2iV ) [ sin ^ J I sin- J (28) 

Again we rewrite this as series in ( 

N-l 



^""^ _|_ ^—m^2^2km _|_ ^(2fc— 2)m _|_ . . . _|_ £ — (2k— 2)m _|_ ^— 2km^^2lm _|_ ^(2£—2)m _|_ . . . _|_ ^— 2£m~j 

71=1 

V-l A; € 

^ ^ ^ ^( 2a + 2 ' 3 + 2 ) m -f 2^( 2a+2 ' 3 ) m _|_ ^(2a+2/3-2)m^j ^9) 



m=l a=— k /3= 



We combine ^( 2 «+ 2 / 3 + 2 ) m i n the first series with 2a 2/3 2 ) m of the last series, and 
C (2o+2/3)m with ^' 2a ~ 2 ^ m . Then (^T]) yields 



a=-k P=-e 



(30) 
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The summand becomes —1 + N(6 a+ p + i i0 + 5 a +/3,o) an d considering separately the 
cases k > £, k < £ and k = £ we find 

2N(2£ + 1) for k > £ 
- (2k + 1)(2£ + 1) + N(Ak + 1) for k = £ (31) 

2N(2k + 1) for k < £ 

This agrees with (|18|). 

5 Higher loop calculations. 

The transition element is now given by 

T(z, y; (3) = (j^yj ' (exp -Ifif*"*) (exp -Is™) (32) 

The factor {g(z)/g(y)} 1 ^ comes from our choice of free and interaction part and our 
normalization of states, in particular, < y\p >= (2 r KTi)~ n l 2 g(y)~ l l 4: exp j-Pjy-' ■ This 
nontrivial measure factor is usually omitted but is crucial to get correct results. 
Vertices are given by 

- ^(3h J\VT + R)dr ; x bg (r) = z + {z - y)r (33) 

and propagators are given in (15-16). One can now compute the loop expansion of 
T; this involves higher loops of a quantum field theory on a finite time segment. Our 
lattice regularization defines all expressions and the results are finite, unambiguous 
and correct (see section 1 for the definition of correct). The two-loop corrections to 
T agree with the result obtained from direct operator methods (see the flow chart). 

6 Phase space path integrals. 

If one moves in the flow chart down on the left hand side, one encounters phase space 
path integrals. Coupling the nN momenta p^j to external sources F\_i (pk lies 
between x fc and and will become equal to zAg fc „!/ 2 /e. We are in the Euclidean 
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case), and the midpoint fluctuations <f k _ x j 2 to Gk-i/2,j one finds, after completing 
squares and integrating out the p's, in the exponent a factor 

E^{-^W + -(?Vi-g£)} 2 (34) 

fc=0 

Expanding this term, one recovers the result already obtained for the discretized 
configureation space path integral, together with an extra F 2 term. It follows that 
the qq and qp propagators in phase space are the same as the qq and iqq propagators 
in configuration space, but the pp propagator is equal to minus the qq propagator 
plus an extra term proportional to 5^/, which cancels the 5^/ present in < qq >. 
Hence, the p propagator is nonsingular. No 5(cr — r) are present in continuum phase 
space Feynman graphs, and the naive approach gives the correct results 



< Pi(a) Pj (r) >= Phgyiz); < g*(<r)^(r) >= -ifiU)(o + 6(r - a)) (35) 

The naive propagator for the kinetic terms | jli[ipiq l — \g l K z )PiPj\d& is given 
( 1 -id a \ ~ 

by G(a,r) = ° 6(a — r). We decompose G as G(a,r) = Gp(a — 

y +id a J 

r)+P(cr, t) where P is annihilated by the field operator (the homogeneous solution) 

and Gp(o — r) — [ 1 ~^ 2 ' ^ ^ ) . The boundary conditions q(a = 

y -^ie(cx - r) A F (a - r) J 

0) = q(a = — 1) = fix P completely, and one recovers (|35| ). Note that one does not 
need any boundary conditions on p, nor is there any need, since there are no zero 
modes in p: all p integrals are convergent and Gaussian. Our discretized approach 
explains this: the variables pk were defined at midpoints (between Xk and Xk-i, and 
were not specified at the endpoints, unlike the Xk for which xn = z and xq = y. 
Similar remarks hold for the ghosts: also they are defined on the midpoints, have no 
boundary conditions and the b, c integrations converge because these are Grassmann 
integrations while the a integrations converge because they are Gaussians. 



7 Mode regularization. 

In this section we will illustrate that the commonly used mode cut-off regulariza- 
tion scheme gives incorrect results for the transition element. In the mode-cut off 
regularization scheme one starts directly from the continuum configuration path in- 
tegral. All quantum fields are expanded in a Fourier series and the path integral is 
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converted into an integral over the Fourier modes. The mode regularization scheme 
amounts to performing all the calculations with a fixed number of Fourier modes, 
say M, and then at the end of the calculation let M — > oo. We shall see that this 
seemingly "natural" regularization scheme is inconsistent with out new Feynman 
rules and therefore yields incorrect result (incorrect in the sense explained in the 
introduction) . 

To be concrete let us consider the same model as before. The continuum action 
is given by 

S = \ /-i l 9ij{Xb9 + q) ^ + > ' ){H ' 1 + ^ + 6V + ^ 

- \pU [ (IT + R)dr ; x bg {r) = z + (z - y)r (36) 

o J —1 

The background fields x\, s (t) satisfy the field equation of the quadratic part = 
■^gij(z) J° 1 (g*g J +&V J +a l a^)d,T and are chosen such that they vanish at the boundary. 

Since the quantum fields q l (r) vanish at the boundary we can expand them in 
the complete set of {sinmrr} on the interval — 1 < r < 0. The ghosts we expand 
into cos niTT since they don't vanish at the boundaries 

oo oo 

Q % = In sin(nvrr); V = ^ b l n cos(nTrr) 

n=l n=0 

oo oo 

c l = c^cos(n7rr); a 1 = a l n cos(n7rr). (37) 

71=0 71=0 

Next we change variables in the path integral from the quantum fields to modes. 
At this stage the measure is fixed by hand such that a Gaussian integral over each 
mode gives one (apart from a possible overall constant). It is straightforward to 
obtain the propagators 

<q l (a)q\r)> = -/3%^A(a,r), (38) 
<b\oy{r)> = -2(3hg ij d 2 a A(a,T), (39) 
<a i {a)a j {r)> = (3hg ij dlA{a,T), (40) 



"° sminna) sinfwrr) 



where 

Note that ( |39l) and (f40|) follow from the identity 2 Y^=i cos niT ' J cos n7! " r + 1 = 
2J2^Li sinnna sinmrr = 5(a—r). (Use that 9{a—r) = — 2J]^Li cosn7ro"sinn7rr/(n7r) 
r is also given by 2 Y^=i sinmrer cosniiT / (mt) + a + 1, and differentiate w.r.t. a). 
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From this identity fl20|) follows. In fact, expanding the ghosts into sines gives the 
same propagators, as the identity shows. 

The propagators < ? l (o")g^(r) > and < q l {(T)q^{r) > and < g l (cr)g J (r) > are 
obtained by simply differentiating fl5B|) appropriately. Mode cut-off regularization 
means that we truncate A (a, r) at some mode M, perform all calculations and at 
the end let M — > oo. 

Let us now illustrate that mode regularization yields incorrect results. Consider 
the two loop graph with contribution 

J = J J drdaA-(a,T) -A(a,r) -A-(a,r), (42) 

where the dot in A(er, r) indicates a time derivative w.r.t. o or r depending on 
which side the dot is (for example A (cr, r) = <9 r A(cr, r)). 

Using ( fi"T| ) and performing the integrals over a and r we get 

1 



1 ' l-(-l) m +"+ fc 



T + 

m + n + k m + n — k 

1 ^ 2 

+ 



m—n+k m—n—k 



(43) 



where the prime indicates that we only sum over m, n, k such that all denominators 
are nonzero. This triple sum is only conditionally convergent. Its result depends 
on the way the summation is performed. Mode cut-off instruct us that we perform 
all sums for a finite upper limit M (the same for all three) and then let the cut-off 
tend to infinity. A numerical calculation yields whereas our Feynman rules give 
Clearly mode cut-off is incorrect for this problem. 



8 Outlook. 

Our results for nonlinear sigma models can serve as a toy model for higher-dimensional 
path-integrals, to clarify there such problems as: equal-time contractions, higher- 
derivative interactions, the measure, boundary conditions, extra ghosts. In partic- 
ular the role of the extra terms due to Weyl ordering is intriguing. If one follows 
Schwinger's analysis of Yang-Mills theory in the Coulomb gauges [10], one is deal- 
ing with a four-dimensional nonlinear sigma model. The operator ordering in the 
Hamiltonian may be fixed by starting with Yang-Mills theory in the Aq = gauge 
(where no ordering ambiguities exist and where it seems therefore reasonable to 
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take the Hamiltonian operator without extra h terms) and then to make a canon- 
ical transformation (at the quantum level!) to the Coulomb gauge. This produces 
extra terms of order h and h 2 in the Coulomb Hamiltonian which Schwinger al- 
ready discovered by requiring that the Poincare generators close at the quantum 
level. According to Christ and Lee [11], Weyl ordering will lead to further % and h 2 
corrections. 

On the other hand, the configuration space approach of Faddeev and Popov 
also ends up with Feynman rules for the same theory, but here there is no sign of 
h corrections and the Feynman rules are straightforward. In fact, the FP approach 
is only intended to yield the Feynman rules at the h — level, but it does not 
address itself to % corrections. Yet, the Coulomb gauge plays a central role in 
fundamental (not practical) discussions of quantum gauge field theory and a precise 
understanding of the quantum theory requires to settle issues at order h and beyond. 
It would therefore be very interesting to generalize our framework and to establish 
a well-defined set of Feynman rules in higher dimensions. The question then arises 
whether our Hamiltonian approach is equivalent to the naive FP method. If not, 
this might have profound implications. 
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